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Abstract. PageRank is a ranking of the webpages used by the search engine 
Google to determine in which order to display webpages. It measures how 
often a given webpage would be visited by a random surfer on the webgraph. It 
seems realistic that this ranking may influence the reputation of the webpages, 
and therefore the behavior of wcbsurfers. We propose a simple model taking 
in account the mutual influence between webranking and wcbsurfing. 

We study this model by considering the following iteration on the set of 
stochastic row vectors: r(s + l) = ui-(r(s)), where for all vector x, u-r(x) is the 
unique invariant measure of a primitive stochastic matrix Mx (x) constructed 
from the adjacency matrix of the webgraph. The parameter T > is fixed and 
measures the confidence of the surfer in the ranking. We call T-PageRank the 
limit, if it exists, of these iterates when s tends to infinity. We also consider 
the simple iteration defined by f (s + 1) = f (s)My (f (s)). 

We prove that, when T is large enough, the fixed point of these iterations 
is unique and the convergence is global on the domain. But for small values of 
T, at least when the matrix C is positive, there are always several fixed points. 

Our analysis uses results of nonlinear Perron-Frobenius theory, Hilbert 
projective metric and Birkhoff 's coefficient of ergodicity. 

1. Introduction 

Motivation. The PageRank algorithm [5] is at the heart of one of the most popular 
Web search engines. Its basic idea is to use the graph structure of the Web in 
order to assign to each webpage a score. Basically, the PageRank score attributed 
to the webpages measures how often a given page would be visited by a random 
walker on the webgraph. Formally, let C = [Cy] be the n x n adjacency matrix of 
the webgraph, so that Cy = 1 if there is a hyperlink from page i to page j and 
Cij — otherwise. For simplicity, we first assume that C is irreducible, that is the 
webgraph is strongly connected. Imagine that, when visiting a page i, a websurfer 
chooses randomly the next webpage he/she will visit, among the pages referenced 
by page i, with the uniform distribution. The trajectory of such a websurfer is a 
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Markov chain with transition matrix M = [M^], given by 



Mi. 



J J2k Cik 

In its most basic version, the PageRank vector r, the entries of which give the 
PageRank score of the webpages, is simply defined as the stationary distribution 
of this random walk on the wcbgraph. Thus r is the invariant measure of the 
irreducible matrix M, that is the unique stochastic vector such that 

r = tM . 

However, the assumption that a websurfer makes uniform draws on the web- 
graph may seem unrealistic: a websurfer could have an a priori idea of the value of 
webpages, therefore favoring pages from reputed sites. The webrank may influence 
the reputation of the websites, and hence it may influence the behavior of the 
websurfers, which ultimately may influence the webrank. In this paper, we propose 
a simple model taking into account the mutual influence between webranking and 
websurfing. 

The T-PageRank. We consider a sequence of stochastic vectors, representing suc- 
cessive webranks, r(0), r(l), . . . , which is defined inductively as follows. The current 
ranking r(s) induces a random walk on the webgraph. We assume that the web- 
surfer moves from page i to page j with probability proportional to Cije E< ^ r ^ j ^ T , 
where E is an increasing function, the energy, and T > is a fixed positive param- 
eter, the temperature. The websurfcr's trajectory is therefore a Markov chain with 
transition matrix M(r(s)), where M(x) is defined for all vector x by 

M(x)-= C - eE(Xj)/T 

The unique stationary distribution of this Markov chain, i.e. the invariant measure 
of the matrix M(r(s)), is then used to update the webrank. Thus, 

(la) r(s + 1) = u(r(s)) , 

where, for all vector x, u(x) is the unique stochastic vector such that 

(lb) u(x) = u(x)M(x) . 

We call T-PageRank the limit of r(s) when s tends to infinity, if it exists, or a fixed 
point of the map u. 

Note that if r(0) is the uniform distribution, then r(l) is the classical PageRank. 
Note also that the temperature T measures the randomness of the process. If T is 
small, with overwhelming probability, the websurfer shall move from page i to one 
of the pages j referenced by page i of best rank, i.e. maximizing r(s)j, whereas if 
T = oo, the websurfer shall draw the next page among the pages j referenced by 
the page i, with the uniform distribution, as in the standard webrank definition. 
For T = oo, the T-PageRank coincides with the classical PageRank, because in this 
case M(r(s)) = M. 

We also consider the simple iteration defined by 

(2) f (s + 1) = f(r(s)) , with f (x) = xM(x) , 

where r(0) is an arbitrary stochastic vector. From a computational point of view, 
this is similar to the standard power method. 
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Main results. Our first main result shows that, if the temperature T is sufficiently 
large, the T-PageRank exists, is unique and does not depend on the initial rank- 
ing. Moreover, at least when the matrix C is primitive, the generalized power 
algorithm ((2|) can be used to compute the T-PageRank. 

Theorem 1.1. If T > nhip(E), where Lip(E') is the Lipschitz constant of the 
function E, then the map u given by \lb\ has a unique fixed point and the iter- 
ates (| la\ converge to it for every initial ranking. Moreover, if C is primitive and 
if T is large enough, the iterates (0) converge to this unique fixed point for every 
initial ranking. 

On the other hand, for small values of T, several T-PageRanks exist, depending 
on the choice of the initial ranking. In some cases, the T-PageRank does nothing but 
validating the initial "belief" in the interest of pages given by the initial ranking. 
Consider for instance the graph given in Figure [T] with adjacency matrix C = 
(J loV and let E(x) = x for all x and T — \. Let r(0) = r(0) = (§ \ +e | - e) 

for an arbitrary small e > 0. Then the iterates |T]) and (2]| converge to a T- 
PageRank close to (0.021 0.978 0.001), so the initial belief that the node 2 is more 
interesting than node 3 has strongly increased. Our second main result shows that 



2 1 3 




Figure 1: For this graph, self-validating effects appear for small temperatures. 



the existence of multiple T-PageRanks is in fact a general feature, when T is small 
enough. 

Theorem 1.2. If C has at least two positive diagonal entries, then multiple T- 
PageRanks exist for T small enough. 

These two theorems follow respectively from Theorems 13.61 and 13. 7\ and from 
Theorem 13. 101 which are stated in a more general framework. 

We also study variants, bringing more realistic models of the behavior of the 
websurfer. This includes the presence of a "damping factor" , as in the standard 
definition of Google's PageRank. We also consider the situation where the websurfer 
may take the webrank into account only when visiting a special page (the search 
engine's webpage). Similar conclusions apply to such variants. 

Method. In order to analyze the map (JTbJ) , we first use Tutte's Matrix Tree Theo- 
rem |37j to express explicitly the invariant measure u(x) in term of the entries of 
M(x). Then we study the convergence and the fixed points of JT]) by using results 
of nonlinear Pcrron-Frobenius theory due to Nussbaum [3T] and Krause [53]. To 
analyze the iteration we use two different approaches, which give convergence 
results under distinct technical assumptions. Our first approach is to show that, 
if T is sufficiently large, the map f from the simplex to itself is a contraction for 
some norm. Our second approach uses Hilbert's projective metric and Birkhoff's 
coefficient of ergodicity. 
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Related work. Several variants of the PageRank are considered in the literature 
in order to have a more realistic model of the behavior of the websurfers. For 
instance, some authors propose to introduce the browser's back button in the 
model [5] [36]. Other try to develop a topic-sensitive ranking by consider- 
ing a surfer whose trajectory on the webgraph depends on his query or book- 
marks [15 ] 120 ] [33] . Moreover, the standard PageRank vector is simply the domi- 
nant left-eigenvector of some stochastic matrix. It has therefore received attention 
from the linear algebra community. Current research concerns for instance efficient 
computation methods [HI [23 HS1 HI 133 , or sensitivity analysis gKHHS] of the 
PageRank. 

Our use of transition probabilities proportional to Cije E ^ s '^' T , where E is 
an energy function and T the temperature, is reminiscent of simulated annealing 
algorithms. For a reference, see Catoni [9]. In the context of opinion formation, 
Holyst et al. [TB] study a social impact model where the probability that an indi- 
vidual changes his opinion depends on a "social temperature" T, which measures 
the randomness of the process. 

The iteration (flj can be studied in the settings of nonlinear Perron-Frobenius 
theory. Many works exist in the literature in order to generalize the classical Perron- 
Frobenius theorems to nonlinear maps on cone satisfying some hypotheses as, for 
instance, primitivity, positivity or homogeneity. See Nussbaum [3T] and also [121125] 
for recent references. 

The iteration f2j has been studied by several authors in an abstract setting. 
Artzrouni and Gavard [3] analyze their dynamics when x(s) behaves asymptotically 
like A s x* for some A ^ 1 and x*. For A — 1, it can be useful to look at the stability of 
a linearization of the system near one of its fixed points [H] . When M(x) is stochastic 
and satisfies certain monotonicity conditions, Conlisk [10] proves the convergence 
of the iterates ([2]) to a stable limit. Lorenz [27] proves their convergence for column- 
stochastic matrices satisfying classic properties of opinion dynamics models. In [28], 
he experimentally studies a reformulation of these models with stochastic matrices 
M(x), where x is an opinion distribution vector. Iterations like could also 
be studied in the setting of nonhomogeneous products of matrices. In this case, 
iterations like x(s + 1) = x(s)Af (s) are considered, where the matrices do not 
depend explicitly on x. Two classical approaches to study their dynamics and 
convergence are the use of ergodicity coefficients (2] [14] [34] [39] or of the joint 
spectral radius [T3] Q31 [15] ■ However, the main results of this paper can not be 
deduced from these works. 

Outline of the paper. We first introduce some preliminaries in Section [2] In Sec- 
tion [3] we prove the main results of this paper, about the existence, uniqueness or 
multiplicity of the T-PageRank, and the convergence of the iterates Jl} and |j2j. 
Then, in Section [4] we introduce a refinement of our model, inspired by the damp- 
ing factor of the classical PageRank algorithm. This variant has the advantage of 
allowing one to work with non strongly connected webgraphs. Moreover, it gives a 
more realistic model of the behavior of a standard websurfer. Section [5] is devoted 
to the estimation of the critical temperature, that is the temperature correspond- 
ing to the loss of the uniqueness of the T-PageRank, for some particular cases. We 
shall see that, even for very small or regular webgraphs, the T-PageRank can have a 
complex behavior. Finally, in Section [6] we experiment the T-PageRank algorithm 
on a large-scale example. 
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2. Preliminaries 

In this section, we introduce some notations and recall some classical concepts 
about nonnegative matrices and graphs, and about iterated maps on cones (for 
background, see [fl E2 EH G2| ) ■ 

We denote by R>o and R>o the sets of nonnegative and positive numbers re- 
spectively. The simplex E = {x 6 R> : ^ Xj = 1} is the set of stochastic row 
vectors. Its relative interior is denoted by int(E) = £ n R> , anc ^ its boundary by 
<9£ = E\int(£). We also write x > y if the vector x — y £ R™ and x > y if the 
vector x — y G R>o- For some vector d, let D = diag(d) be the diagonal matrix 
such that Du = for all i and = for all i ^ j. Finally, let 1 denote the 
column vector of all ones, and Sij denote the Kronecker delta. 

Let us begin with nonnegative matrices and graphs [HEIGH]. Let A £ R™q™ 
be a nonnegative matrix. Its associated directed graph G(A) is a graph with nodes 
1,2, ... ,n, and a directed edge if and only if A t j > 0. This graph is strongly 
connected, that is there exists a directed path between each pair of nodes, if and 
only if A is irreducible. It is moreover aperiodic if and only if A is primitive. Let 
r be a node of G(A). A directed subgraph R of G(A) which contains no directed 
cycles and such that, for each node i ^ r, there is exactly one edge leaving i in R, 
is called a spanning arborescence rooted at r. The set of spanning arborescences of 
G(A) rooted at r is denoted by A(r). 

The transition probabilities of a Unite Markov chain, or equivalently a random 
walk on a directed graph with weighted edges, can be represented by a stochastic 
matrix, that is a matrix A £ R™q" such that Al = 1. A stochastic vector u £ E is 
an invariant measure of A if uA = u. The Perron-Frobenius theory ensures that 
an irreducible stochastic matrix has a unique invariant measure, which is moreover 
positive. An expression of this invariant measure is given by Tutte's Matrix Tree 
Theorem [37] (for a proof, see Section 9.6 of [7J). 

Theorem 2.1 (Matrix Tree Theorem, Tutte). Let A £ R"p™ be an irreducible 
stochastic matrix, and let u be its invariant measure. Then u = v/J^Vj, where 
for all index r 

(3) Vr = Yi n a « ■ 

ReA(r) (i-J)eR 

We now present some useful concepts about iterated maps on cones [311 E] • A 
function h is subhomogeneous on a set U if Ah(x) < h(Ax) for all x £ U and every 
< A < 1. It is order-preserving on U if h(x) < h(y) for all x, y £ U such that 
x < y. Hilbert's projective metric du, is defined as 

d H ■ M> x M" -> M>o : (x, y) i-> max In . 

This metric defines a distance on int(E). An induced projective metric between 
two positive matrices can be defined as 

d H : K™q™ x M™q™ — M> : (A, B) ^ sup d H {xA, xB) . 
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The coefficient of ergodicity tb, also known as Birkhoff's contraction coefficient, is 
defined for a nonnegative matrix A having no zero column as 

fA . d H (xA,yA) 

TB (A) = SUp 



<Mx,y) 

x#Ay ' 



x,y6K" 



This coefficient tb(A) G [0, 1], and tb(A) < 1 if and only if A is positive. 

The following theorem deals with the fixed points and the convergence of iter- 
ated maps on a cone. This theorem is a simple formulation of the very general 
Theorems 2.5 and 2.7 (or more precisely of Corollaries 2.2 and 2.5) of Nussbaum 
in [31]. 

Theorem 2.2 (Nussbaum). Let h: K" — > R™ be a continuous, order-preserving 
map which is subhomogeneous on int(E). Suppose moreover that if x. G int(E) is 
an eigenvector of h, then h is continuously differentiable on an open neighborhood 
of x and the matrix h'(x) is nonnegative and irreducible. Then h has at most 
one eigenvector in int(E). As a consequence, the map u: int(E) — > int(E): x 
h(x) /y~lj h(x), has at most one fixed point. 

Moreover, ifh has an eigenvector x G int(E), and i/h'(x) is primitive, then all the 
orbits of u converge to its unique fixed point x. 

We end with this preliminary section with a result of Krause [23j . also dealing 
with the fixed point and the convergence of iterated maps on the positive orthant. 

Theorem 2.3 (Krause). Let p: R> — > R>o be a continuous map which is not 
identically and such that p(Ax) = Ap(x) for all x G R> , A > 0; and p(x) < p(y) 
for all 0<x<y. Let U = {x G R£ : p(x) = 1}. 

Let h: R™ — > R™ be a map such that there exist a, (3 > and v G R™ such that 
av < h(x) < /3v /or a// x G U. Suppose also that Ah(x) < h(y) for all x,y 6 {/ 
and < A < 1 such that Ax < y; and that Ah(x) < h(y) if moreover A < 1 
and Ax ^ y. Then h has a unique eigenvector in U . As a consequence, the map 
u: U — * int(L7): x h(x)/p(h(x)) has a unique fixed point in U. Moreover, all 
the orbits of u converge to this unique fixed point. 

3. Existence, uniqueness, and approximation of the T-PageRank 

3.1. Hypotheses. Let C be an n x n irreducible nonnegative matrix. For all 
temperature T > 0, and all x G E, let My(x) be the irreducible stochastic matrix 
such that 



M T (x) 



where gr ■ [0,1] — > R>o is a continuously differentiable map with <7t(0) = 1. We 
suppose moreover that gx is increasing with g^ : [0, 1] — > R>o and make the follow- 
ing assumptions on the asymptotic behavior of gx 

(A ) lim L mhl (g T ) = oo , 

(Aoo) lim L max (g T ) = , 

T — >oo 

where L m in(ffT) and imax(S'T) are defined as 

/ n • 9t( x ) i t , s 5t( 2: ) 

W 5r = mm — TT and ^ma X (5T) = max — ^ . 

as o,i] ffTl^J *€[o,i] sr(a;) 
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Example. These hypotheses are satisfied in particular for g(x) = e E ( x " T , where 
the energy function E is continuously differentiable, with E(0) — and E'(x) > 
for all x G [0,1]. In this case, L min (g T ) = mm xe[0 E'(x)/T and L max (g T ) = 
maxagp.i] E'(x)/T. 

3.2. Preliminary results. The following elementary lemmas will be useful in the 
sequel. 

Lemma 3.1. For all x,y G [0,1], 

Lminigr) (x - y) + < (in gT ^M < L max (g T ) (x - y) + , 
V 9T{y)J 

where for all x + = max{0,a;}. Moreover, if x,y =/= 0, 

ln^) + <W 5 T) fln^ + . 

9T{y)J V vJ 

Proof. Let x,y G [0,1]. Then 

(In g T (x) - In g T (y)) + > min 9t ^ (x - y) + = L mia (g T ) (x - y) + , 
ae[o,i] g T (a) 

Qng T (x) - In g T (y)) + < max ilJL ( x _ y )+ = L raax ( 3T ) (x - y) + . 
oe[Q,i] Jr(a) 

Moreover, (x — y) + < (In x — In j/) + if x, y G ]0, 1] . □ 

Lemma 3.2. Let x, y £ [0, 1]. Ifx>y, then lim^^o 5t(^)/5t(j/) = °°- 

Proof. This result follows directly from Lemma 13.11 and Assumption flAoH ■ D 



Lemma 3.3. TTie map (?t tends to the constant function equal to 1, uniformly in 
[0, 1], w/ien T tends to infinity. 

Proof. For all a; G [0, 1], by Lemma |3~T1 

grp (x\ 

lng T (x) = In — — < L max (g T ) (x - 0) < L max (g T ) ■ 
9t(0) 

Therefore, by Assumption flA^p | , 

lim sup higrix) < lim LmaxCffr) = . □ 

T ^°°2;G[0,l] T ^°° 

It follows from Lemma [3.31 that limT->oo Afr(x) = diag(Cl) _1 C uniformly for 
x G S. 

3.3. Fixed points and convergence of u^. When C is irreducible, for all x G S, 
the matrix Myfx) is irreducible, and we can define the map 

ut : S — » S : x i— > (x) , 

that sends x to the unique invariant measure Uy(x) of Mr(x). The Matrix Tree 
Theorem enables us to express explicitly ut(x). 
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Lemma 3.4. Assume that C is irreducible. Then, the invariant measure o/Mr(x) 
is given by 

h T (x) 



u T (x) 



E fe h T (x) fe ' 



wh 



ere 



(4) h T (x) r =(^C r fe5T(x fe ))( Yl IT C »9T^j) 

k R.eA(r) (i,j)eli. 

Proof. Apply Theorcm l2.1l to M^(x), and take h^(x) = [i v, where v is given by ([3]) 
and (j. = UiJ2k C ik9T(^k)- □ 

The existence of fixed points for is then proved using Brouwer's Fixed Point 
Theorem. 

Proposition 3.5. Assume that C is irreducible. The map ut has at least one fixed 
point in int(E). Moreover, every fixed point of Ut is in int(S). 



Proof. By Lemma l3.4[ the map u^: E — > E is continuous, and therefore Brouwer's 
Fixed Point Theorem ensures the existence of at least one fixed point for uy. 
Moreover, since the invariant measure of an irreducible matrix is positive, and 
since Mr(x) is irreducible, ut sends E to int(E), and therefore every fixed point 
of ut is in int(E). □ 

The following result concerns the uniqueness of the fixed point and the conver- 
gence of the orbits of u T . With Assumption jA^l satisfied, it shows that the map 
Uj- has a unique fixed point and that all its orbits converge to this fixed point, for a 
sufficiently large temperature T. It can be proved using Nussbaum's Theorem l2.2[ 
as we do it here, or, under the same hypotheses, using Krause's Theorem 12.31 (take 
P( x ) = E,: x 0- 

Theorem 3.6. Assume that C is irreducible. If nL max {gj') < 1, the map has 
a unique fixed point x^, which belongs to int(E). Moreover all the orbits of Ux 
converge to this fixed point. 

Proof. Since g-r is increasing, h<r is an order-preserving map from R" to itself: 
x < y implies hj^x) < hr(y). Now, let us show that is subhomogeneous on 
int(E). Let x G int(E) and < A < 1. Any entry of hr(x) is a sum of positively 
weighted terms like 



k 

with ^ fe 7fe = n. By Lemma 1331 for each k G {1, . . . , n}, 

In — tt r < L max (g T ) In - . 

5T(Ax fc ) A 

Therefore, if nimaxCffr) < 1, then A 1 '™^^) < gr(^k), and it follows that 
Ahvr(x) < hvr(Ax). Since < hvr(O), this shows that h*r is subhomogeneous on 
int(E). 
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Finally, let xt be a fixed point of ut, by Proposition l3.5l The derivative h T (x) 
is a nonnegative continuous function of x: 



dh T (x) 



C r i g' T (xg 



ReA(r) (i,j)eR 



k ReA(r) 9T(xe) ( . j}( _ R 

where m^i? = (i,^) G Moreover, since g T takes positive values, dh.T(x) r /dxi > 
as soon as C r e > or me t R > for some R G A(r). In particular, vn r ^ > for all 
R G ■A(r). Since an irreducible matrix with positive diagonal is primitive, h T (x) is 
also primitive (see for instance Corollary 2.2.28 in [4]). Therefore, by Theorem l2.2[ 
xt is the unique fixed point of ut, and all the orbits of ut converges to Xt- □ 

3.4. Fixed points and convergence of fy. We now consider the map 

f T : E -> E: x i-> xM T (x) . 

The following result shows that if T is sufficiently large, the fixed point of the map 
ut can be computed by iterating fy. 

Theorem 3.7. The fixed points o/u^ and fx are the same. If C is primitive, then, 
for T sufficiently large, all the orbits of fp converge to the fixed point xt of u-r . 

Proof. Clearly, fy and uy have the same fixed points. Suppose now that C is 
primitive, and let us show then that, for T is sufficiently large, fy is a contraction 
for some particular norm. For every x, y G E and for any norm || • ||, 

||f T (x)-f T (y)|| <sup||(x-y)f T (v)|| . 
The derivative of fr satisfies 



— - M TMtl + X (% - ^<v), j; StCa!)T(vt) 
By Assumption flApd ) and Lemma [ 



lim f T (v) = lim M T (v) = diag(Cl)~ 1 C , 

uniformly for v G E. Let A = diag(Cl)" 1 C, and let S={zel": J^k z k = 0} be 
the space of vectors orthogonal to the vector 1. The map x 1— > xA preserves the 
space <S, because Al = 1. Let As denote the restriction of the map x 1— > xA to S. 
Since the matrix A is primitive, its spectral radius, p(A), is a simple eigenvalue and 
all the other eigenvalues of A have a strictly smaller modulus. Moreover, the space 
S contains no eigenvector of A for the eigenvalue p(A), because such an eigenvector 
u must be a scalar multiple of the Perron eigenvector of A, which has positive 
entries, contradicting J2 k u^. = 0. We deduce that p(As) < p(A) = 1. It follows 
that there exists a norm || • || such that |||As||| < 1, where ||| • ||| is the matrix norm 
induced by || • || (see for instance Lemma 5.6.10 in [T7]). Therefore, since f T (v) 
tends uniformly to A for v G E when T tends to 00, 



lim S up|||f T (v) s |||HIM<l 
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where fy(v)s denotes the restriction of fp(v) on S. It follows that, for all a G 
] HI As 1, 1[, there exists T a such that for all T > T Q , 

||f T (x)-f T (y)|| < ||x-y|| sup|||f^(v) 5 ||| <a||x-y|| . 

Hence, for such temperature T, by Banach's Fixed Point Theorem, fy has a unique 
fixed point and every orbit of f-r converges to this fixed point. □ 

Example. For the convergence of the orbits of fr, the primitivity of C cannot be 
dispensed with. Let for instance C = { \ J) and T > 0. Then Mf(x) = C for all 
x G E, and fy(x) oscillates when k tends to infinity, unless x = (i A). 

For positive matrices C G let us now derive another convergence criterion, 

depending on Birkhoff 's coefficient of ergodicity. 

Lemma 3.8. Assume that C is positive. Then, for any x, y G int(E), 
d H (M T (x),M T (y)) < 2L max (g T )d H (x,y) . 

Proof. Let x, y G int(E) be fixed. Let us define a = maxjj ^Tpp 2 an d = 



mmi J' Mr(y)lj - By dennition ' 

tl/f(Mr(x),Mr(y)) = sup maxln 



(zM T (x)). t (zM r (y))j 
«eK» D "*j"" (zAf T (y)) 2 (zM T (x))j 



< sup maxln ^^7^% =lng 
- zG R5 « (zM T (y)), 09zM T (y))j /3 

Moreover, by Lemma |3.1[ 

Cij 9t fa ) E fc Cjfc g T (y fc ) 



In a = max In 



/ 1 5t(xj) 5r(yfc) 
< max m — + In 



i,fe V 9r(yj) ' 5T(x fe ) 

< L max (g T ) max ( In ] + max [ In — 
\ i \ y 3 J k V x fe 

= L max (g T ) ( [ max In — j + ( max In — 

\V yjJ \ k x fe 

= L max (g T )d H (x,y) , 

since x,y G int(E) implies maxj In 2i > and maxfe In — > 0. We get similarly 
- In (3 < L max (g T ) d H (x, y) . □ 

Proposition 3.9. Assume that C is positive. If 2L max (gr) < 1 — tb(C), then 
has a unique fixed point G int(E) and all the orbits of fr converge to this fixed 
point. 



Proof. Let x, y G int(E). By Lemma [37 

dff(f T (x),f T (y)) < 4(xM T (x), y M T (x)) +rf H (yM T (x),yM T (y)) 

< r B (M T (x))d H (x,y) + d ff (M T (x), Af T (y)) 

< (t b (<7) + 2i max ( 5T )) d H (x,y) • 
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Therefore, if 2L max {gT) < 1 — tb{C), then fr is a contraction on int(S) with respect 
to the distance dn- By Banach's Fixed Point Theorem, fr has a unique fixed point 
xj- £ int(S) and all the orbits of fr converge to this fixed point. □ 

3.5. Existence of multiple fixed points of uj- and fy. Theorems 13.61 and 13.71 

show that for a sufficiently large temperature T, the maps and fr have a unique 
fixed point. We can naturally wonder about the uniqueness of the fixed point of 
these maps for small T: we show that, at least when C is positive, multiple fixed 
points always exist. 

Theorem 3.10. Assume that C is irreducible and that the first column of C is 
positive. Then, for all < e < ^, there exists T e such that for T < T e , the map ut 
has a fixed point in S E = {x £ S, xi > 1 — e}. 

Assume now that C is irreducible with C\\ > only, and that there exists e n > 0, 
independent of T , such that 5T(£n)"~ 1 < 5t(1 — £n) for all T > 0. Then, for all 
< e < e n , there exists T £ such that for T < T e , the map ur has a fixed point in 
X £ = {xe S,xi > 1-e}. 

Proof. Let k be the number of indices i ^ 1 such that Cn > 0. Let < e < \, 
and let x G S e . By the irreducibility of C, there exists a spanning arborescence R 
rooted at 1, containing all the k arcs (i, 1) with i / 1 and Cn > 0. Hence 

h T (x)i > Cn5 T (x 1 ) J] CijgTixj) > ag T {^i) k+1 > ag T (xi) fc ffT (l - e) , 

where a = C\\ \\u ^ eR dj > 0. Let r ^ 1. If C r \ ^ 0, then a spanning arbores- 
cence rooted at r can have at most k — 1 arcs («, 1) with Cn > 0, whereas it can 
have at most k arcs («, 1) with Cn > in general. Hence, in all cases, hr(x) r is 
a sum of positively weighted terms like WigT^t) 11 with ^2 e "fe — n and 71 < k. 
This implies that, for r =/= 1, 

h T (x) r < (3g T ( Xl ) k g T (e) n - k , 

for some positive constant (3. Therefore, 

. . 1 1 
u ^Wi = — — r-7=r- > 



1+V h T (x),. ~ 1 , / in/? g T (e)"- k 
1 + l^r^l h T (x)i + ^J Q g T (l- e ) 

If the first column of C is positive, then k = n — 1. By Lemma 13. 11 

In g^ 1 - £ ) > L min ( gT ) (1 - 2e) . 
grie) 

Therefore, if L m i n (5r)(l - 2e) > l n (" M we ge t u T (x)i > 1 - e. This shows 
that ut(S £ ) C S e . By Brouwer's Fixed Point Theorem, the continuous map uy 
has therefore at least one fixed point in E e . 

Now, suppose we know only that C\\ > 0, but there exists e n > such that 
5t(£ti)" _1 < 5t(1— £n) for all T > 0. The map : e 1— > (n— 1) ln^yfe)— ln<7;r(l — e) 
is increasing and its derivative satisfies <p' T {e) > nL m i n (<rr)- Let < e < e n . We 
have ifT{sn) — Vt(^) > fi£min(<7T)(£n — £)■ Moreover, fc > 1 by irreducibility of C 
and </?T(£n) < 0, hence 

. pr(l — e) . Or(l — e) , N , , w ., 

5t(£)™ * 3t(£)' 1 1 
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Therefore, if nL min {gT){e n — e) > In — we get ut(x)i > 1 — e. The result 

follows by the same argument as above. □ 

Remark 3.1. When <?t(ie) = eP^ x >l T for some increasing energy E, then e n satis- 
fies the condition gT^n ) n_1 < g T (l-e„) for allT>0 if and only if (n-l)E(e n ) < 
E(l — e n ), which holds for some < e n < 1, since E(0) — and E(l) > 0. 

Corollary 3.11. If C is positive, then, for T > sufficiently small, the map ut 
has several fixed points in E. 

Example. If C is not positive, the existence of several fixed points for small T is not 
insured. Indeed, we shall see in Remark |5. II that for C = {\ g) and gr{x) — e x l T , 
the fixed point of and fy is unique for each T > 0. 



4. Refinement of the model 

In the present section, we study a more general model, which includes a damping 
factor < 7 < 1, as is the standard definition of Google's PageRank : 6\, in which 
the websurfer either jumps to the search engine with probability 1 — 7 or moves 
to a neighbor page with probability 7. The presence of a damping factor yields a 
more realistic model of the websurfer walk. Moreover, it allows one to deal with 
reducible matrices, and it improves the convergence speed of iterative methods. 

Let C be a n x 11 nonnegative matrix with no zero row and let d £ K™ be a 
personalization vector. For all temperature < T < 00, let us define as previously 
<7t: [0, 1] — > M>o as a continuously differentiable and increasing map, with <7t(0) = 
1 and g' T : [0,1] — ► M>o- Suppose that Assumptions JAqD and flApd ) are satisfied. 
For a temperature T = 00, let us also define g^ (x) = 1 for all x G [0, 1]. 

For any two temperatures < Ti,T^ < 00, and for all x 6 E, we can consider 
the positive transition matrix Mt 1j t 2i7 (x) defined as 

l^k C ik 9Ti (x fe ) 2^fc d fc 9t 2 (x fe j 

Remark 4.1. for simplicity, we consider the same family of weight functions gr 
for the first and the second term of Mr lj j' 2!7 (x). Note however that the results of 
this section remain true if two families gr ± and gT 2 are considered. 

Remark 4.2. Suppose T\ = 00, Ti < 00 and < 7 < 1, and let x be the current 
ranking vector. Then Mj- 1: t 2 ^ 7 (x) is the transition matrix of the following random 
walk on the graph. At each step of his walk, either, with probability 7, the websurfer 
draws the next page uniformly among the pages referenced by his current page. Or, 
with probability 1— 7, he refers to the Web search engine, and therefore preferentially 
chooses for the next page a webpage with a good ranking. 

The maps UTi.Tj/y an d £11,22,7 are defined as previously: ut , 1 ,t 2 ,7(x) is the unique 
invariant measure of Mt^Ts^x) and fT lt T 2 ,-y(x) — xAfT 1 ,T 2 ,7(x). 

Theorem 13.61 about the uniqueness of the fixed point of ut can be adapted in 
the following way. 

Proposition 4.1. If nL max (g Tl ) + (n - l)L max (g T2 ) < 1, the map UTi,Ta,7 has a 
unique fixed point x.r ly T 2 ,-y i n E. Moreover, all the orbits o/u^ j^ ^ converge to the 
fixed point x^^^. 
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d ^f )r =c«9' T Mi){ e n 

1 ReA(r) (i,j)eR 

^*.w)( E ( n E 



Proof. If Ti = T2 = 00, the result follows directly from Perron-Frobenius theory. 
Let us therefore suppose that Ti < 00 or T2 < 00. For every x G E, by Theorem l2.1| 
u Ti,t 2 , 7 = h Tl! T 2i7 (x)/X; fc h T ^ T2 ^(x) fe , where 

h Tl>T2 , 7 (x) r = (^C rfcffTl ( Xfc ))( ^ [] W( X )i 

with PF(x).y = J2 k (7 C-ij 5Ti (Xj) d fc g T2 (x fe ) + (1 - 7) C jfc ff Tl (x fe ) ffT2 (x 3 )) . 

Since and gr 2 are nondecreasing, h^ ^ ^ is an order-preserving map. More- 
over, assume that nT m ax(<?Ti) + (n — l)imax(ffTi) < 1- Then, as in the proof of 
Theorem 13.61 hT 1 ,T 2 ,7( x ) is shown to be subhomogeneous on int(E). 

Finally, the derivative h' Ti T ^ 7 (x) is a nonnegative continuous function: 



k K ReA(r) (i,j)£R (i,j)£R 

where 

dW Q^ 13 = S V E ( 7 Cu g Ti ( x ^) dfc 9t 2 (x fe ) + (1 - 7) C ik g Tl (x fc ) d £ 3 ^ 2 (x^)) 
+ 7 Cy .97i ( x j) d £ (x £ ) + (1 - 7) C u g' Ti (x e ) dj gr 2 (xy) . 

Let us now prove that Ta 7 ( x ) is a positive matrix for every x. Suppose first 
that Ti < 00. Then g' T2 (x £ ) > 0, and there exists a spanning arborescence R £ A(r) 
and a node i such that (i,^) £ i?, since G(Mj 1i t 2 , 7 (x)) is the complete graph. It 
follows that, for this R, 

and hence 9hT 1 ,r 2 ,7( x )r/^ x f > 0- Now, suppose that T2 = 00 and Ti < 00. Then 
we can suppose without loss of generality that C has no zero column (see Remark l4.3l 
below). Either C r i > 0, and therefore dh^ (x) r /9x £ > 0. Or there exists i ^ r 
such that Cm > 0, and for all R £ A(r), there exists j such that £ R, that is 

^^>(l-7)C«d i ^ I (x,)>0 , 

and hence 9hT 1 ,T 2 ,7( x )r/9x £ > 0. 

Since Brouwer's Fixed Point Theorem ensures the existence of at least one fixed 
point XTi,T 2 ,7 £ int(E) for the continuous map u^,^^ which sends E to int(E), 
by Theorem 12.21 this fixed point x^^^ is the unique fixed point of 117^, t 2 , 7; and 
all the orbits of u^^^ converge to xt^t^- D 



Remark 4.3. If T2 = 00 and the matrix C has a zero column, the problem can 
be reduced to a problem of smaller dimension with a matrix with no zero column. 
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Indeed, suppose the nth column of C is zero. Then 

, i •> r \ \ Efc d fe - (1 -7) d n - 

(, u Ti,oo, 7 lXji ••• Ut-^qo 7 (3C) n _iJ — 117^,00, 7W ) 

Z^fc d fe 

UTi,oo,7(x)t, = (1 - 7) ' 

where x = (xi ■•• x„_i) and Ut 1j00i7 (x) is £/ie invariant measure of a matrix 
-^Ti,oo,7(x), loii/i C £/ie principal submatrix of C corresponding to the indices 
1, . . . , n — 1, and d some positive vector of length n — 1. 

The following adaptations of Theorem 13 . 71 and Proposition ^. 91 about the unique- 
ness of the fixed point of fj> are quite direct. 

Proposition 4.2. The fixed points of ^Ut 1 ,t 2 ,7 flndfri,T2, 7 are the same. Moreover, 
for Ti and T 2 sufficiently large, all the orbits 0/^^2,7 converge to the fixed point 

X Ti,T 2 , 7 0/u Tl ,T 2 ,7- 

Proposition 4.3. Assume that C is positive. If 2(T max (<7^ 1 ) + I ma x(te)) < 
1 — rs(C), i/ien fri,T 2 ,7 a unigue fixed point xt 1 ,t 2 , 7 G int(S) and aZZ i/ie orbits 
o/iTi,T 2 ,7 converge to the fixed point t 2 ,7- 

We next show that the map fri,T 2 ,7 has multiple fixed points if either Ti or T2 
is sufficiently small. Of course, this will require the damping factor to give enough 
weight to the terms corresponding to the small temperature in equation (|5j). 

Proposition 4.4. Assume that the first column ofC is positive and that i < 7 < 1. 
Then, there exists Tq > such that, for all T± < Tq and for all Ti S [0,oo], the 
map fTi,T 2 ,7 has a fixed point in {x G S : xi > i}. 

The conclusion of Proposition 14.41 is weaker than that of Theorem 13.101 The 
latter shows that for a sufficiently small temperature, we can find a fixed point of 
fr arbitrarily close to a vertex of the simplex, whereas the former shows that for 
a sufficiently small temperature, we can find a fixed point of fri,T 2 ,7 m the region 
Xi > 1/2 of the simplex. In fact, such a fixed point may not approach a vertex of 
the simplex as the temperature tends to 0. This discrepancy is due to presence of 
the damping factor. The proof of Proposition. 14.41 is a straightforward adaptation 
of the proof of Theorem 13.101 

Remark 4.4. If the first column of C is not positive, the existence of a fixed 
point such that xi > x 2 , . . . , x„ for small T\ is not insured. Consider for instance 
C = (11), 9t{x) — e x l T , T2 = oo, d = (1 1) and < 7 < 1. Then, for each each 
T\ > 0, any fixed point xt 1j0 o,7 o/fT lj0 o,7 belongs to [0, |[ x ]i, 1]. 

We now show the existence of multiple fixed points for a sufficiently small tem- 
perature T 2 . 

Proposition 4.5. Let < e < \ and let /1 = dfc/di. If0<j<e, there 

exists T £i7 such that for T2 < T e-7 and for all T\ € [0, 00], the map foo,T 2 ,7 has a 
fixed point in S E = {x e S : xj — 
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Proof. Let x e E e . Since by Lemma I3TT] ln(<?y 2 (l — £)/.9T 2 (e)) > £imn(<?T 2 ) (1 — 2e), 
we have 

fT T ,( X ) 1 = 7 V ^ni^j +d- 7 ) d ^( x i) >(i_ 7 ) dl ^( Xl ) 

1-7 



> 



„ e -i mto ( 9 T 2 )(l-2e) 



Therefore, if L min (g T2 )(l-2e) > In , then f Tl ,T 2 , 7 (x)i > l-£ and f 00! T 2i7 (x) G 

S s . By Brouwer's Fixed Point Theorem, the continuous map fri,T 2 ,7 has at least 
one fixed point in S E . □ 

We now consider the case where the damping factor 7 approaches 1 and T\ = 
00. Then, the corresponding value of the generalized PageRank converges to an 
invariant measure of the matrix diag(Cl) _1 C, independently of the choice of T2. In 
order to prove this, we will need the following classical result, which is a particular 
case of Corollary 3.1 in [3D]. For a matrix M G C nx ™, the index, indM, is the 
smallest nonnegative integer k such that rankM + = rankM fc , and its Drazin 
inverse, M D , is the unique solution X of the equations M k+1 X = M k , XMX = X, 
MX = XM, where k = indM. 

Lemma 4.6 (See [3D], Coro. 3.1). Let M G C nxn . IfbxdM = 1, then\im E ^ e(M+ 
el)- 1 =1 -MM D . 

Proposition 4.7. Assume T\ — 00. For all vector norm \\ ■ \\ and all e > 0, there 
exists 7 e < 1 such that for every fixed point x 0/ foo t 2 -y, with j e < 7 < 1 7 there 
exists an invariant measure u of diag(Cl) _1 C such that ||x — u|| < e. 

Proof. Let A = diag(Cl) _1 C, and let M = I— A. Since A is stochastic, indM = 1 
(see for instance Theorem 8.4.2 in [3]). Therefore, with 7 = (1 + e) -1 , 

limfl -7) (J - 1A)- 1 = lim e(M + eiy 1 = I - MM D . 

Let || • || be a vector norm and ||| ■ ||| its induced matrix norm, and let v > such 
that ||v|| < v for all stochastic vector v. Let e > 0. There exists 7 e < 1 such that 
if 7 £ < 7 < 1, 

|||(1 - 7 ) (J - 1 A)- 1 -{I- MM D )\\ < v- x e . 

Let 7 G ]7 e , 1[, and let x be a fixed point of foo,T 2 ,7j that is x = v(x)(l — "f){I — 
jAy 1 , where v(x)< = d { g T2 (x<)/X) fc d fe gr 2 (xfc) for all i. Then, 

||x -v(x) (7 -MM D ) || < e . 

But v(x)(7 — MM D ) is an invariant measure of the matrix A. Indeed, I — MM D 
is stochastic, and (I - MM D )(I - A) = M - MM D M = M - M 2 M D = 0, by 
definition of the Drazin inverse. □ 



5. Estimating the critical temperature 

We call critical temperature the larger temperature for which the number of 
fixed points of ut changes. It corresponds to the loss of the uniqueness of the fixed 
point. In this section, we are interested in estimating the critical temperature for 
some particular cases. We will study in details the case of n x n matrices of all ones 
with the particular weight function gr{x) = e x ^ T . 
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We suppose that gr(x) = e x l T and first consider the particular case where the 
graph is complete with 

/l ••• V 



c = 



g 



v ••• h 

For this matrix, the point x = i 1 is a fixed point of for all T. We are interested 
in the existence of other fixed points, depending on the temperature T. 

Lemma 5.1. Assume that C is the n x n matrix of all ones. The point x is a fixed 
point of ut if and only if x £ E and there exists A £ R such that A = Xie _Xi / T /or 
every i = 1, . . . , n. 

Proof. This follows directly from x = fr(x). □ 

Lemma 5.2. Assume £/ia£ C is the n x n matrix of all ones. The point x is a fixed 
point of ut if and only if x 6 E and there exists K, C {1, . . . , n}, y £ [0, T] and 
z >T such that ye~ y l T — ze~ z l T , \K\y + (n — |/C|)z = 1, Xj = y for all i £ K, and 
Xj = z /or aZ/ i 1C. 

Proof. Since the map x i— > xe~ x / T is increasing for < a; < T and decreasing for 
x > T, there can be at most two values y z such that ye~ y l T — ze~ z l T — A for 
a given A G M. The result hence follows from Lemma I57T1 and x £ E. □ 

Since in our case, L max (gT) = T^ 1 , we know from Theorem 13.61 that the critical 
temperature is at most n. Proposition 15.31 shows that this critical temperature is 
in fact roughly (Inn)" 1 when n tends to infinity. 

Proposition 5.3. Assume that C is the n x n matrix of all ones. If n > 2, the 
map ut has a unique fixed point x =(—•••—) j/ and only ifT> T*{n), where 

1 - -3- 1 - i i 

— < T*(ri) = sup — f < 



ln((lnn- l)n+ 1) ~ v ' a> i ln((a - l)n + 1) ln(n-l) ' 

thus T*(n) ~ when n tends to oo. If n = 2, the map Ut has a unique fixed 
point x = (i i) if and only if T > T*(2) = 

Proof. From Lemma T5.21 x G int(E) is a fixed point of Uj, with (— • ■ • -), if 
and only if there exists K, C {1, . . . , n}, y, z £ R, such that Xj = y for i G /C, = z 
for i ^ /C, < k = \K\ < n, ky+(n — k)z — 1, ye~ y l T — ze~ z l T ' , and y < z. Denote 
a = Since y < i, we get necessary that a > 1. From ye~ y l T = ze _z / T , we get 
7" = ITa.fci where 

1-i 



(»-*)i»(^i 2 +0 



This implies that ut has a fixed point x G int(E), x ^ ( ) ^ an( i on ly if 

T G T = {r aifc ,a > 1,A; G {l,...,n- 1}}. Let 

1-i 

T (n) = supT Q ,„_i = sup— a . 

q>i a>i ln((a - l)n + 1) 

We shall show that T = ]0, T*(n)} when n > 2 and T = ]0, T*(2)[ when n = 2. 

First, a study of T a ^ as a function of k shows that it is increasing. It is therefore 
sufficient to show that {T a)7l _i, a > 1} = ]0,T*(n)] when n > 2, and {T ai i,a > 
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1} = ]0, T*(2)[ when n = 2. Second, a study of T Q)Tl _i as a function of a > 1 shows 
that, when n > 2, T Q ,„_i is increasing, then decreasing, tends to when a goes 
to infinity, and its maximum is attained for a = a n , where a n > — — — . Hence 
T = ]0, T*(n)]. When n = 2, T Qj i is decreasing, tends to when a goes to infinity, 
and to \ when a goes to 1. Hence T*(2) = i, and T = ]0, T*(2)[. 
Moreover, for n > 3, 

1- J- 1 

T*(n) = T Qni „_i = — - -f — — — < 



In((a„ - l)n + 1) ln((a„ - l)n + 1) ' 

and since a n > 2 ^ 1 ^ , we get T*(n) < ln („ 1 „ 1 ^ ■ For the lower bound, we get 
21nn,n-i < T*{n), since Inn > 1. □ 



Proposition 15.31 deals with the very special case of a complete graph. In more 
general circumstances, the exact computation of the critical temperature seems 
out of range. However, we can obain numerically a lower bound of the critical 
temperature, which seems to be an accurate estimate, using the following homotopy- 
type method. We first choose two random initial vectors on the simplex. Then, we 
iterate the map fy from each of these vectors. For small values of T, this yields 
with an overwhelming probability two different webranks. Then, we increase the 
temperature T, and keep iterating the map fr on each of these webranks, until 
the two webranks coincide. This yields a lower bound of the critical temperature. 
Then, we repeat this procedure, with new random initial vectors, until the lower 
bound of the critical temperature is not improved any more. Note that the simpler 
method consisting in keeping T fixed and iterating fy from various initial conditions 
(random vectors or Dirac distributions on a vertex of the simplex) experimentally 
yields an under estimate of the critical temperature. 

Using the previously described homotopy-type method, we computed numeri- 
cally the critical temperature for two families of graphs. These experiments reveal 
that the l/ln(n) asymptotics obtained for the complete graph gives a good gen- 
eral estimate. We first considered the ring graph, with n nodes, in which node 
i is connected to its two neighbors and to itself. The critical temperature, for 
n = 51,201,501 and 1001 is shown by stars in Figure O The exact value of the 
critical temperature of the complete graph with n nodes, T*(n), is drawn as a 
continuous curve. We see that the critical temperatures of the ring and complete 
graphs are essentially proportional. We also computed numerically the critical 
temperature for a standard model of random directed graph, in which the presence 
of the different arcs are independent random variables, and for every the 
probability of presence of the arc is given by the same number p. We took 

p = 10/n, so that every node is connected to an average number of 10 nodes. The 
corresponding critical temperatures are represented by circles. The values of these 
critical temperatures do not seem to change significantly with the realization of 
the random graph, hence, each of the values which are represented correspond to a 
unique realization. 

Remark 5.1. Let us briefly discuss about the case of an arbitrary 2x2 irreducible 
matrix C with the weight function gx = e x ^ T . In this case, some elementary cal- 
culations give information about the critical temperature pQ. Firstly, the critical 
temperature for a graph of only two nodes is always less than 1, since it can be 
shown that uj- has a unigue fixed point if T > 1 . This is the best general upper 
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Figure 2: Estimation of the critical temperature T c as a function of the number of nodes n for 
the complete graph (continuous curve), the ring graph (stars) and random graphs (circles). 



p. a/10 



0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

T T 
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Figure 3: For 2x2 matrices, the map has 5 fixed points for (T, a) in the black region, 3 fixed 
points for (T, a) in the grey region, and 1 fixed point otherwise. 



bound that can be given for problems of this dimension since for every T < 1, we 
can construct a 2 x 2 matrix such that ut has at least two fixed points. Moreover, 
one can show that, for every T > 0, the map ut has at most 5 fixed points and 
does not have any orbit of period greater than 1. Numerical experiments show that 
for a 2 x 2 irreducible matrix, the number of fixed points of the map ut can change 
0, 1, 2 or even 3 times when decreasing the temperature T. 

This can be seen on Figures which were obtained experimentally. Let a — 
C11/C12 and (3 = C22/C21. For a specified (3, vlt has 5 fixed point if (T, a) belongs 
to the black region, 3 fixed points in the grey region and 1 fixed point in the white 
region. 

6. Experiments on a subgraph of the Web 

In this section, we briefly present our experiments of the T-PageRank on a large- 
scale example. We consider a subgraph of the Web with about 280000 nodes which 
has been obtained by S. Kamvar from a crawl on the Stanford wetQ. We use the 
variant of our model presented in Section |4j with a transition matrix given by 

C--e^l T e x ^ T 
M(*k= 7 ^ % ^ + (1-7) 



oXfc /T 



where we suppose that for each dangling node i (i.e. a node corresponding to a 
webpage without hyperlink), the ith row of the matrix C is a row of all ones. 



^The adjacency matrix can be found on http: //www. Stanford. edu/~sdkamvar/research. html. 
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The chosen damping factor is 7 = 0.85. Wc have computed the T-PageRank from 
the recurrence ([2]) for various temperatures T and initial rankings. As expected, 
when the temperature T is large, the T-PageRank is very close to the classical 
PageRank, and when T approaches zero, arbitrary close initial rankings can induce 
totally different T-PageRanks. The critical temperature experimentally seems to 
be about T — 0.033. It has the same order of magnitude as the T*(n) = 0.06148 
estimate discussed in Section [5l 

As in [38] , we represented in a log- log scale the cumulative distribution function 
of the PageRank, i.e the proportion of pages for which the T-PageRank is larger 
than a given value, as a function of this value. In Figure (a), we show the suc- 
cessive T-PageRanks obtained for increasing temperatures from T = 0.015 to the 
critical temperature T = 0.033 by the following variant of the previously described 
homotopy method: for T = 0.015, we iterate the map f-r, with a Dirac mass on a 
vertex of the simplex as initial ranking, until a fixed point is reached. Then, for 
each new value of T, we iterate fr until a fixed point is reached, starting from the 
previous fixed point. For T < 0.032 the distribution of the T-PageRank is quite 
different from this of the PageRank and it comes closer suddenly for T = 0.033. In 
Figure (b), we show the successive T-PageRanks obtained by a similar method for 
decreasing temperatures from T = 0.033 to 0.009, with the classical PageRank as 
an initial ranking. The latter procedure may be compared with simulated annealing 
schemes, in which the temperature is gradually decreased. Until T = 0.0091, the 
distribution of the T-PageRank is quite similar to that of the classical PageRank 
(see a zoom in Figure (c)). With T = 0.009, the T-PageRank moves suddenly 
away from the PageRank. These figures suggest that the gap between webpages 
considered as "good" and "bad" is more pronounced with the T-PageRank than 
with the classical PageRank. 

We have also compared the five best nodes for the classical PageRank and for 
the T-PageRank with decreasing temperatures T = 0.033, 0.015 and 0.0091. As 
we see in Table [TJ for T = 0.033, the PageRank and T-PageRank give a similar 
ranking for the top-five. But for smaller temperatures as T = 0.015 or T = 0.0091, 
even the two best nodes are exchanged. 

T = 0.033 T= 0.015 T = 0.0091 



1 2 2 

2 1 1 

3 3 3 

4 6 46 

5 7 33 



Table 1: The five best nodes of the T-PageRank for several values ofT: the numbers refer to the 
rankings acccording to the classical PageRank. 

Since for this special set of data, the correspondence between the page numbers 
and the urls is not available, one cannot interpret the discrepancies between the 
PageRank and the T-Pagerank in Table 6.1. In [32 , J. -P. Poveda made similar 
experiments on the larger matrix obtained by S. Kamvar for a crawl of the union 
of the Stanford and Berkeley Webs3, with about 685000 nodes, for which, this 



'Also available on http://www.stanford.edu/~sdkamvar/research.html. 
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Figure 4- Fraction of pages having a PageRank larger than a particular value, (a) T-PageRanks 
computed with increasing temperatures from T = 0.015 to 0.033, with a vertex of the simplex 
as initial ranking, (b) T-PageRanks computed with decreasing temperatures from T = 0.033 to 
0.009, with the classical PageRank vector as initial ranking, (c) Zoom of Figure (b). 



time, the correspondence between some pages and the main urls is given. These 
experiments suggest that the T-PageRank obtained by the latter scheme, in which 
the temperature is gradually decreased, as illustrated in Figure (b), might be of 
practical interest. A full experimentation, on a real scale web, is beyond the scope 
of the present paper. 

As a concluding remark, we would like to point out that our results might be 
considered as an argument in favor of the thesis that one should not use PageRank 
type measures to assess quality. Indeed, the validity of the classical PageRank 
relies on an ideal view of the web, in which the makers of pages are thought of as 
experts, creating hyperlinks only to pages they carefully examined, and judged by 
themselves to be of interest. In the real world, however, the web makers may be 
influenced by factors like reputation, to which the webrank participates. Within 
the limits of the model, our results show that pathological phenomena, like getting 
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non unique and even not meaningful rankings, may occur in a self-referential world 
in which the websurfers would excessively rely on the web ranking, rather than on 
their own judgement. 
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